The center-of-mass response of confined systems 

L. F. Lemmens 

Dept. Natuurkunde, Universiteit Antwerpen RUCA, Groenenborgerlaan 171, B2020 

Antwerpen, Belgium. 

F. Brosens and J. T. Devreese 
Dept. Natuurkunde, Universiteit Antwerpen UIA, Universiteitsplein 1, B2610 

Antwerpen, Belgium. 
(September 11, 1998, resubmitted December 22, 1998 to Phys. Rev. A) 

For confined systems of identical particles, either bosons or fermions, we argue that the parabolic 
nature of the confinement potential is a prerequisite for the non-dissipative character of the center 
of mass response to a uniform probe. For an excitation in a parabolic confining potential, the 
half width of the density response function depends nevertheless quantitatively on properties of the 
internal degrees of freedom, as is illustrated here for an ideal confined gas of identical particles with 
harmonic interparticle interactions. 
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I. INTRODUCTION 

In the present paper we study the response of a confined system of spin-polarized identical particles if a uniform but 
time-dependent force is applied. Our method is based on the Feynman-Kac functional representing the propagator 
of distinguishable particles, followed by the projection of this functional on the irreducible representations of the 
permutation group as required for identical particles. The projection introduced in |5J] is applied on confined systems. 
Details have been given earlier P| to derive the thermodynamical properties and the static response functions for 
bosons and for fermions. An application of the method to real systems can be found in Ref. || for rubidium. 

The response function is calculated for a spatially homogeneous but time-dependent force. It is shown that if the 
confinement potential is a polynomial of order two, the motion of the center of mass does not induce transitions between 
the modes representing the relative motion of the internal degrees of freedom. This follows from the independence of 
these degrees of freedom, a property known || for harmonic systems. The analogue in extended systems is the Kohn 
theorem Q , well known for the cyclotron resonance in an interacting electron gas. The connection to the response of 
center-of-mass excitation becomes especially relevant if the two-body interaction only depends on the distance vector 
between the positions. If one of these conditions is not satisfied, we shall show that the homogeneous external force 
couples to the internal degrees of freedom. Even if both conditions are satisfied the shape of the response function still 
depends on the characteristics of the internal degrees of freedom, typically the frequency of the internal oscillation 
modes and the number of particles in the well. This suggests that these characteristics can be estimated from the 
shape of the response to a homogeneous time-dependent force. For instance, monitoring the density of a confined 
system whilst changing the direction of the gravitational force on the center of mass or changing slightly the form of 
the confinement potential gives information on the internal degrees of freedom and the number of particles confined 
in the system. 

In view of the experimental progress made ||-^| on the Bose-Einstein condensed systems, the present response 
property of a perfect harmonic confined system could be relevant to check the parabolicity of the trap. Further recent 
experiments fllCfl use strong deviations from the parabolic confinement potential to probe the time dependent behavior 
of the confined atoms in a condensed phase. In this respect it seems important to know the implications of perfect 
parabolicity. Of course, the response of the system to such a local perturbation of the density would be of major 
importance. 

The paper is organized as follows. In Sec. II the formal response theory is developed. In Sec. Ill we illustrate the 
response properties for the exactly soluble model of a confined gas with harmonic interparticle interactions. In the 
last section a discussion and the conclusions are given. 
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II. PARABOLICITY AND THE KOHN THEOREM 



In a 3 AT dimensional configuration space with Yj (j — 1,2, ... , N) denoting the positions of the N particles, the 
position of the center of mass is defined as R = r j- For a quadratic confinement potential the center of mass 

can be introduced as follows: 



JV N 



(1) 



From mathematical statistics it is well known that for Gaussian distributions the average R can take any value 



without affecting the value of J2jLi ( r j ~ ^-Y: i- e - t ne center of mass could be independent of the deviations from that 
center. This remains true in a physical system when the two-body potential depends only on the distance between 

particles, because Ylj i—i V ( r j ~~ r i) — Ylj i=i V (( r i ~ R<) ~ ( r 2 ~ R<)) • I n the harmonic model the following two-body 
potential is used: 
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This two-body potential allows for exact solutions of the propagator and the projection on the symmetric or antisym- 
metric irreducible representation of the permutation group can be carried out. 
Introducing a homogeneous time-dependent force: 



JV 



f(r).^r J= iVf(r).R 



(3) 



the potential confining N particles and the interaction of the particles with each other and the external field are 
specifying the system for which the response will be calculated. Atomic units with h — m — 1 are used throughout 
this paper. 

If the particles were distinguishable, the Euclidean-time propagator Kjj (r", rlr', 0) for this interacting system can 
be obtained from a straightforward generalization of our calculation in Ref. [0. Denoting by r S M. 3N the set of 
position vectors ri, . . . , rjv, the result turns out to be: 



K d (T",t\T',0) = 



A^(VATR",T|ViVR',0 



\/jVf(7 



3=1 



(4) 



K (VNR",t\VNR',0 

where K (r r , r|ro, 0)1^ denotes the well-known propagator of a three-dimensional harmonic oscillator of frequency w: 

v (rl + i-q) cosh 



tf(r T ,T|r o ,0)l 



27r sinh i 



exp 



2r T • r 



sinh 'i 
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The quantity K (R", r|R, 0)[^ (t)] denotes the propagator of a three-dimensional harmonic oscillator of frequency 
fi in the presence of a driving force f (r) . If the driving force acts in the real time t, the complex time r = j3 + it has 
to be introduced, with (3 = 1/kgT with kg the Boltzmann constant and T the temperature. This propagator can be 
obtained in a direct way using the functional integration techniques illustrated in B leading to the following result: 



A-(R",r|R,0)| f o (T) =K(K",r\R',0)\ i 



(6) 
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f (s) • f (a) sin (fl (t — s)) sin (ft (s — a)) dsda x 
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ds 



^From the structure of the propagator Ko (r", r|r', 0) for distinguishable particles, it is clear that projection on the 
representations of the permutation group will not affect the quotient of propagators that contains the center of mass. 
Therefore, the propagator Ki (r", r|r', 0) for identical particles is given by 
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Kj (f",r\f',0) 



K (//VR",t|%/]VR / ,0 



li (r",r|r',0), 



(7) 



where Kj (r", r|r', 0) accounts for the permutations of the particles, and is discussed extensively in |2j. Although the 
present derivation has been given for the fully harmonic model (harmonic confinement and harmonic interaction), 
the line of the derivation remains valid for more general two-body interactions. Indeed, the separation between the 
propagation of the center of mass and of the other degrees of freedom is unaffected by the introduction of the two-body 
potential if it only depends on the difference vectors r j — r; . 



III. RESPONSE TO AN UNIFORM FORCE 



Knowing the propagator of the interacting many-particle system, the density response can be calculated in a 
straightforward way as the following expectation value: 

1 N ( 1 \ 3 f 

= jvS<*( r - r i)) = ( ^ J / <(*)ex P Hq-r)dq, 

i=i 
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qW = ]v2Z( eX P( iC l- r i))- 



(8) 
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The averages (A (r, t)) are defined as: 

{A (r, *)) = J dvA (r) Kj (r, ft + it\ffi) , 

where Z is the partition function for the system without external force at inverse temperature p. Performing this 
averages means the repetitive calculation of Gaussian integrals leading to the following result: 



(9) 



(t) =<(0)exp -iq- / f(s) 



n f (0) = exp 
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w sinhiw/3 risinhifi/3 

The response function n f (r, t) is then obtained by a Fourier transform leading to a Gaussian density distribution 
function 
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(11) 



in which the variance not only depends on the frequency f2 of the confining potential, but also on the frequency w of 
the internal degrees of freedom: 



o2 1 [N- 1 , wfi 
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(12) 



In the absence of interparticle interactions, one readily obtains S 2 \ w= n — (cothifi/?) / (2J1) , which in the low- 
temperature limit ft — > oo gives the variance of an harmonic oscillator with frequency fl. 
In the average density response 

(13) 



(r (t)) = J rn f (r,t)dr = -^J^ f (s) sin (0 (t - s)) ds 



one clearly sees the resonant structure in the presence of an oscillating driving force f (s) = f v sinus with a specific 
frequency v: 

f v Q sin tv ~ v sin fit 



(r (*)) |f( s)= f oi 



n 



Q 2 -v 2 



(14) 



This expectation value of the position does not depend on the parameters of the internal degrees of freedom, as 
required by the Kohn theorem. 
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IV. DISCUSSION 



From the derivation of the density response to a uniform time-dependent force, it is clear that the factorization (Q) 
is crucial to obtain a linear relation between the density and the applied force that does not depend on the internal 
degrees of freedom. Higher-order terms in the confinement potential will spoil this relation. A term of the third degree 
in the confinement potential already introduces a contribution V con f ~ 3R J^Lj (ij) 2 , that mixes the center-of-mass 
motion with that of the internal degrees of freedom, and makes transitions possible in the internal degrees of freedom 
leading eventually to dissipative behavior of the center-of-mass excitation. 

Furthermore the Gaussian nature of the density response may allow one to check the parabolic character of the 
well by measuring the mean density and comparing it with the trap frequencies. The standard deviations of these 
measurements contain information about the excitation frequencies of the internal degrees of freedom and the number 
of particles contained in the well. The response of a center-of-mass mode behaving according to the Kohn theorem 
contains information about the internal degrees of freedom. Its resonance structure reveals the parabolic quality of 
the well. 

It should be noted that the statistics of the internal degrees of freedom does not enter in the expressions. Therefore 
the response of a parabolic well to a uniform external force does not distinguish between fermions or bosons merely 
because it only applies tot the center of mass. 



V. ACKNOWLEDGMENTS 



This work is performed in the framework of the FWO projects No. 1.5.729.94, 1.5.545.98, G. 0287.95, G.0071.98 and 
WO.073.94N (Wetenschappelijke Onderzoeksgemeenschap over "Laagdimensionele systemen" ) , the "Interuniversitairc 
Attractiepolen - Belgische Staat, Diensten van de Eerste Minister - Wetenschappelijke, Technische en Culturele 
Aangelegenheden" , and in the framework of the BOF NOI 1997 projects of the Universiteit Antwerpen. One of the 
authors ( F.B.) acknowledges the FWO (Fonds voor Wetenschappelijk Onderzoek-Vlaanderen) for financial support. 



[1] R.P. Feynman, Statistical Mechanics A set of Lectures, W. A. Benjamin Reading, 1972 

[2] F. Brosens, J. T. Devreese, and L. F. Lemmens, Phys. Rev. E 55, 227 (1997); E 55, 6795 (1997); E 57, 3871 (1998); E 
58, 1634 (1998). 

[3] J. Tempere, F. Brosens, L. F. Lemmens, and J. T. Devreese, Solid State Commun. 107, 51 (1998). 
[4] W. Kohn, Phys. Rev. 123, 1242 (1962) 

[5] A. Skyrme and R.J. Elliott, Proc. Roy. Soc. A 232, 566, (1955) 

[6] M. H. Anderson, J. R. Ensher, M. R. Matthews, C. E. Wieman, and E. A. Cornell, Science 269, 198 (1995). 
[7] K. B. Davis, M.-O. Mewes, M. R. Andrews, N. J. Van Druten, D. S. Durfee, D. M. Kurn, and W. Ketterle, Phys. Rev. 
Lett. 75 , 3969 (1995). 

[8] C. C. Bradley, C. A. Sackett, J. J. Tollett, and R. G. Hulet, Phys. Rev. Lett. 75, 1687 (1995). 
[9] C.C. Bradley, CA. Sacket, and R.G. Hulet, Phys. Rev. Lett. 78, 985 (1997). 
[10] D. M. Stamper-Kurn, M. R. Andrews, A. P. Chikkatur, S. Inouye, H.-J. Miesner, J. Stenger, and W. Ketterle , Phys. Rev. 
Lett. 80, 2027(1998) 



4 



